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ABSTRACT 

This paper provides a fast and patient-specific scatter artifact correction method for cone-beam computed tomog¬ 
raphy (CBCT) used in image-guided interventional procedures. Due to increased irradiated volume of interest 
in CBCT imaging, scatter radiation has increased dramatically compared to 2D imaging, leading to a degra¬ 
dation of image quality. In this study, we propose a scatter artifact correction strategy using an analytical 
convolution-based model whose free parameters are estimated using a rough estimation of scatter profiles from 
the acquired cone-beam projections. It was evaluated using Monte Carlo simulations with both monochromatic 
and polychromatic X-ray sources. The results demonstrated that the proposed method significantly reduced the 
scatter-induced shading artifacts and recovered CT numbers. 
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1. INTRODUCTION 

Due to increased irradiated image volumes, scatter radiation has been one of the major challenges in fiat- 
panel detector based cone-beam computed tomography (CBCT). Typical visual appearances of scatter artifacts 
include cupping, dark streaks between dense objects, degraded contrast resolution, and inaccurate CT numbers 
for quantitative measurements. Many scatter correction methods have been proposed since the early days of 
x-ray imaging. These include scatter rejection technique^ such as the use of an air gap, antiscatter grid and bow- 
tie filter, analytical modeling of scatter profiles ,1^®^ Monte Carlo simulations,primary modulation methods 
using either a spatial primary modulatoiElllini or a temporal primary modulatoi!^ and other measurement-based 
scatter correction methods which measure scatter radiatiorP^M^ using either a beam stop,l^ an aperture,!^ or a 

While these methods have been demonstrated to mitigate scatter artifacts, some challenges still remain. 
For example, scatter rejection techniques usually yield inadequate correction and additional scatter corrections 
are needed.l^ Scatter rejection techniques can also reduce soft-tissue contrast-to-noise ratio (CNR) due to the 
concomitant attenuation of primary radiations.!^ In current analytical modeling method, both calibration mea¬ 
surements and Monte Carlo (MC) simulations may be needed to calibrate the free parameters in the analytical 
modeling of scatter profiles. MC simulations can be computationally expensive and may take prohibitively long 
processing time. In primary modulation methods, the necessary extra hardware support or mechanical modifica¬ 
tions add complexities to the current image acquisition system and thus become difficult to adapt into a clinical 
setting. In direct scatter measurement methods, extra scans or extra hardware or a dedicated reconstruction 
algorithm may be needed. As a result, it is highly desirable to develop a scatter artifacts correction method that 
is fast, patient specific, sufficiently accurate and, ideally, does not require additional scans, or extra hardware 
support and modifications. In this paper, we present a new method to achieve this goal. 


2. METHODS 

The flowchart of the proposed algorithm is presented in Figure 1. The proposed algorithm estimates scatter 
radiation distribution using a coarse scatter distribution which is regularized with a convolution-based scatter 
model. The coarse scatter estimation is achieved by taking the difference between the raw projection data and the 
polychromatic reprojection data of a segmented image volume. The method includes three major components: 





Figure 1. Flowchart of the proposed scatter correction method. 

convolution-based scatter estimation model, coarse scatter estimation and calibration of free parameters used in 
the convolution-based model. 

In the convolution-based scatter model, scatter signal can be expressed as a convolution between a 
scatter potential ^ and a scatter kernel 


I^^{a,c,d) = $[ p ( q :), co , ci ] * K{di,d2), 


( 1 ) 


where ^[p((a),c] = cq + ci • and K{d) = ^-di{x-d2)‘^ with c = (co,ci) and d = (^1,^2) 

respectively. Here p denotes the raw projection and x denotes the detector index. The constant cq represents the 
contributions from the Compton scattering while the second term was attributed to the Rayleigh scattering.!^ 
Although the scatter radiation distribution can be estimated using the above convolution-based model effectively, 
a calibration experiment is often needed to determine the free parameters sets c, d. Methods have been proposed 
to determine these free parameters using either Monte Carlo simulations or a dedicated calibration scan. In this 
paper, we employed a novel coarse scatter estimation method to determine these free parameters. 

The raw projection data I was modeled as a summation of primary beam Ip and scatter contribution A, 


I — Ip Is 


( 2 ) 


Instead of directly estimate A, we generated an estimate of the primary signal Ip which was estimated to be the 
forward projection of segmented image volume using a polychromatic x-ray spectrum. To estimate the scatter 
contribution, we simply subtracted the estimated primary signal Ip from the measured raw projection data I to 
obtain the needed coarse estimation of the scatter signal Ig. 

The estimated scatter profile A was used in the analytical convolution model in Eq. Q to calibrate the 
unknown parameters c and d. By minimizing the least square error between Ig and I^^ ^ the unknown parameters 
could be obtained. Using the obtained final analytical convolution model, the final scatter radiation distribution 
was calculated and the corrected projection data was then obtained by subtracting the final scatter signal from 

































Figure 2. Scatter correction for the quasi-thorax phantom with 70 keV monochromatic X-ray source. After correction, 
scatter artifacts are almost completely removed and CT number are comparable to the primary image. The difference 
images show each image minus the scatter-free primary image. Display window: [-1200HU, 500HU] for both the CT 
images and the difference images. 


the raw projection data. Using the scatter corrected projection data, the routine FBP reconstruction method 
was used to obtain CT images with reduced scatter induced artifacts. 

3. NUMERICAL SIMULATIONS 

In this study, an anthropomorphic quasi-thorax phantonQwas used to generate both primary beam and scatters 
using a Geant4 based Monte Carlo simulation package, GATE.^^Our numerical simulation studies contained two 
aspects: in the first numerical simulation study, a monochromatic x-ray beam at 70 keV was used since beam 
hardening artifacts have similar profiles as scatter artifacts. To be realistic, a polychromatic x-ray spectrum at 
100 kVp was used in our second numerical simulation study. These two studies were evaluated to demonstrate 
the proposed method, since the polychromatic forward projection model can be used to take the beam hardening 
contribution in account during the coarse scatter estimation. The fan angle and the cone angle of the simulated 
data acquisitions are 23.6° and 16°, respectively. The geometrical configuration of the simulated GBCT system 
was similar to the Varian 2100 EX System used in image-guided radiation therapy. Scatter projection, primary 
projection, and total projection were recorded separately. During the correction procedure, the Monte Garlo 
simulation toolkit Geant4 was applied to generate the polychromatic X-ray spectrum for the 100 kVp case. 

4. RESULTS 

4.1. Monochromatic case 

Fig-i shows the results of a scatter-free primary image, an image with both primary and scatter contributions, 
and the corrected GT image using the proposed algorithm. To help visualize the performance of the algorithm, 
the difference between the aforementioned three images and the scatter-free primary image was calculated and 
shown in the second row in Eig. . In the corrected image, scatter artifacts are accurately removed and the 
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Figure 3. Line profiles of the quasi-thorax phantom with monochromatic 70 keV X-ray source. Shading and reduced 
CT number caused by scatter radiation are clearly shown by the red line. After correction with the proposed method, 
shading artifacts are well removed and CT number are well recovered, and the corrected prohles match the scatter-free 
primary prohles well, (a) Vertical direction, (b) horizontal direction. 

CT number is restored to be closer to the primary image. The average errors in CT number were reduced from 
-190 HU to -11 HU for the different ROIs (shown in Figure. upper-left corner). Fig. shows the line prohles 
of the quasi-thorax phantom along central vertical direction and central horizontal direction (shown in Figure.]^ 
upper-left corner). Scatter caused shading artifacts are almost completely removed and the reduced CT number 
are successfully recovered. The corrected prohles match the scatter-free primary prohles very closely. 

4.2. Polychromatic case 

Fig.i shows the results of the same phantom but with the 100 kVp polychromatic x-ray spectrum. As shown 
in the corrected image and the corresponding difference image, scatter artifacts were also signihcantly reduced 
and the CT number is comparable to the primary image. The average CT number errors are reduced from -206 
HU to -23 HU for the different ROIs (shown in Figure. upper-left corner). Fig. shows the line prohles of 
the quasi-thorax phantom along central vertical direction and central horizontal direction (shown in Figure. 
upper-left corner). We can see scatter caused shading artifacts are almost completely removed and the reduced 
CT number are accurately recovered. 

5. CONCLUSION AND FUTURE WORK 

This work presents a new approach to correct scatter artifacts for the CBCT imaging systems. In the proposed 
method, instead of using a dedicated calibration measurement or Monte Carlo simulations to calibrate the free 
parameters in the convolution model, we used the difference of the raw projection data and the polychromatic 
reprojection as our estimated rough scatter prohle to calibrate the free parameters in the analytical convolution 
kernel. The proposed method was evaluated using Monte Carlo simulations with both monochromatic and 
polychromatic X-ray sources. Results demonstrated that scatter artifacts were almost completely removed in 
both cases. Significant improvements have been achieved in both image uniformity and CT number accuracy. 
Note that the proposed correction algorithm is intrinsically patient-specific. The proposed method avoids the 
need of extra hardware support, system modifications, or prolonged data acquisition time. Thus the proposed 
method can be potentially employed to correct scatter artifacts for existing CBCT scanners. 

There are several analytical scatter estimation models. For example, scatter signal can be expressed as 
weighted summation of a set of sine and cosine basis function. The proposed method can be applied to these 
analytical scatter estimation models to calibrate the unknown weight. Further research will be focused on other 
scatter models. 




















Figure 4. Scatter correction for the quasi-thorax phantom with 100 kVp polychromatic X-ray source. After correction, 
scatter artifacts are almost completely removed and CT number are comparable to the primary image. The difference 
images show each image minus the scatter-free primary image. Display window: [-1200HU, 500HU] for both the CT 
images and the difference images. 



Figure 5. Line prohles of the quasi-thorax phantom with polychromatic 100 kVp X-ray source. Shading and reduced 
CT number caused by scatter radiation are clearly shown by the red line. After correction with the proposed method, 
shading artifacts are well removed and CT number are well recovered, and the corrected profiles match the scatter-free 
primary prohles well, (a) Vertical direction, (b) horizontal direction. 
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